First-order nonthermal phase transition after preheating 
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During preheating after inflation, parametric resonance rapidly generates very large fluctuations of 
scalar fields. In models where the inflaton field (f> oscillates in a double-well potential and interacts 
with another scalar field X, fluctuations of X can keep the <f> — * symmetry temporarily restored. 
If the coupling of (f> to X is much stronger than the inflaton self-coupling, the subsequent sym- 
metry breaking is a first-order phase transition. We demonstrate the existence of this nonthermal 
phase transition with lattice simulations of the full nonlinear dynamics of the interacting fields. In 
particular, we observe nucleation of an expanding bubble. 
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Cosmological phase transitions are one of the central 
topics of modern cosmology |lj . Recently, this theory was 
supplemented by the possibility of nonthermal cosmolog- 
ical phase transitions i.e. phase transitions driven by 
fluctuations produced so rapidly that they did not have 
time to thermalize. Large nonthermal fluctuations natu- 
rally occur in inflationary models during preheating ||. 

Fluctuations of Bose fields generated by the paramet- 
ric resonance during preheating have large occupation 
numbers and can be considered as interacting classical 
waves, which allows to study the dynamics of fluctua- 
tions during and after preheating using lattice simula- 
tions Numerical calculations, as well as analytical 
estimates ]4[-|ll| have shown that the maximal values 
achieved by fluctuations can be large enough to cause 
cosmologically interesting phase transitions. Neverthe- 
less, until now there was no direct demonstration of the 
existence of such phase transitions; several groups which 
studied this issue numerically [ fl2[ either concentrated 
on models where phase transitions cannot occur, or ne- 
glected essential feedback effects such as rescattering of 
created particles. 

We have performed a number of lattice simulations of 
nonthermal phase transitions, which demonstrated for- 
mation of various types of topological defects Jl3|]. In 
this Letter we report results that prove that nonthermal 
phase transitions may take place after preheating even 
on the scale as large as the GUT scale ~ 10 16 GeV. The 
phase transition that we have found is first-order, which 
may have particularly important cosmological implica- 
tions. First-order phase transitions have a very clear 
signature: they proceed through nucleation and subse- 
quent expansion of a bubble of the new phase inside the 
old phase. In our opinion, the presence of this distinc- 
tive signature eliminates all doubts about the possibility 
of nonthermal phase transitions in the class of theories 
under investigation. 

As a prototype we will use the model 
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The inflaton scalar field <j> has a double- well potential and 
interacts with an A^-component scalar field X; X 2 — 
Y^Li Xf. For simplicity, the field X is taken massless 
and without self-interaction. The fields couple minimally 
to gravity in a FRW universe with a scale factor a(t). 

The initial conditions at the beginning of preheating 
are determined by the preceding stage of inflation. One 
can define the moment when preheating starts as the 
moment when the velocity of the field (f> is zero in con- 
formal time rj, drj = dt/a(t), a(0) = 1. This happens 
when (f)(0) w 0.35Afpi ||. We will use rescaled confor- 
mal time r = y/\(f>(0)ri. Inhomogeneous modes of <f> and 
all modes of X are taken to be in their conformal vacua 
at r = 0. We consider the case v < 10~ 3 Afpi, g 2 S> A, 
and take A = 10 -13 The strength of the resonance 
depends nonmonotonically on the resonance parameter 
q = ,9 2 /4A, being maximal around q = n 2 /2 |10|. The 
condition <7 2 /A ~ 2n 2 ^> 1 means that the evolution be- 
gins in the regime of a broad parametric resonance. 

The equations of motion for (f> and X are 

4> + 3H(f> - V 2 q)/a 2 + AO 2 - v 2 )(f) + g 2 X 2 (f, = , (2) 
Xi + 3HXi - V 2 X t /a 2 + g 2 tfX l = , (3) 

where H = a/ a. Substituting (f> — (fo + 8(f), where (f>o is 
the inflaton's homogeneous ("zero") mode, into (0) and 
assembling terms linear in <fro, we obtain the effective 
mass of d>: 
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Angular brackets denote spatial averaging. Because fluc- 
tuations rapidly become classical [Q, their variances can 
be computed as spatial averages; thus (X 2 ) is the vari- 
ance of X. The maximal value of (X 2 ) grows with TV, 
while that of ((5(f)) 2 ) = (((f) - (f) ) 2 ) does not §,|l|. Wc 
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therefore expect that at large N fluctuations of X will 
play a more important role in @j than those of <fi, at 
least for some time. In general, there is no guarantee 
that a useful effective potential can be defined for a state 
far from thermal equilibrium. We have found that in 
our case Eq. (Q) gives a good working definition for the 
effective mass. 

The idea of nonthermal phase transitions [g] is that 
large fluctuations of (X 2 ) (and ((S(j)) 2 )) generated during 
preheating can change the shape of the effective poten- 
tial and lead to symmetry restoration. Afterwards, the 
universe expands, the values of (X 2 ) and ({5(f>) 2 ) drop 
down, and the phase transition with symmetry breaking 
occurs. In the case g 2 /X ^> 1 the phase transition can 
be of the first order. This can be expected by analogy 
with the usual thermal case jli]]. Let us establish the 
necessary conditions for this transition to occur and to 
be of the first order. 

(i) At the time of the phase transition, the point <p — 
should be a local minimum of the effective potential. 
From (||), we see that this means that g 2 (X 2 ) > Xv 2 . 

(ii) At the same time, the typical momentum of X 
particles should be smaller than gv. This is the condition 
of the existence of a potential barrier. Particles with 
momenta p < gv cannot penetrate the state with \<f>\ « v, 
so they cannot change the shape of the effective potential 
at 1 <fi\ w v. Therefore, if both conditions (i) and (ii) are 
satisfied, the effective potential has a local minimum at 
<p = and two degenerate minima at <p w ±u. 

(iii) Before the minima at <f> ss ±v become deeper than 
the minimum at <f> = 0, the inflaton's zero mode should 
decay significantly, so that it performs small oscillations 
near </> = 0. Then, after the minimum at \4>\ « v becomes 
deeper than the minimum at <f> = 0, fluctuations of <fi 
drive the system over the potential barrier, creating an 
expanding bubble. 

Before performing numerical investigation of any par- 
ticular model one may want to find out in which cases 
these conditions can be satisfied. If the initial value of </>o 
is much larger than v, the first stages of preheating pro- 
ceed in the same way as in the conformal theories (v = 0), 
which is a well established case. The maximal value of 
(X 2 ) is reached at some time r max . The previous stud- 
ies §,|ll| indicated that (X 2 ) max ~ v^ 2 (0)/ga 2 (r roax ). 
We need, however, an estimate for (X 2 ) at the time of 
the phase transition. If, at r > r max , (X 2 ) merely de- 
creased with time due to the redshift, we could estimate 
it by simply replacing a(r max ) with the current value of 
a. In reality, fluctuations of X decrease somewhat slower, 
because of the continuing decay of the inflaton, but since 
the condition (i) holds whenever v is small enough com- 
pared to fluctuations, it is safe to make the replacement 
a(r max ) — > a. As for the condition (ii), we note that the 
scale of momenta of X particles is set by the frequency of 
the inflaton's oscillations, p» ~ V~X(f>(0) / a, although nu- 
merical factors are to be expected in this estimate. Both 
conditions (i) and (ii) are met if ^P- < v 2 < VN^-. 



We see that the window allowed for v 2 grows with q and 
with N. Therefore we will explore cases g 2 /X'^> 1 and 
N>1. 

The condition (iii) cannot be verified using only results 
of the previous studies. Nevertheless, we may expect a 
fairly rapid decay of the zero mode by analogy with such a 
decay in the model with a massive inflaton B . The anal- 
ogy is relevant when the amplitude of (f>o in the present 
model becomes comparable to v, so that deviations from 
the conformal invariance set in. This expectation is well 
confirmed by our numerical results. 

For numerical studies, the full nonlinear equations of 
motion (||), (||) were solved directly in the configuration 
space. The computations were done on 64 3 lattices for 
N= 1, 2, 9 and on 128 3 lattices for N = 1, 2. Below we 
present results for the model with parameters g 2 / X = 200 
and d = 0.7x 10" 3 M P1 w 0.8 x 10 16 GeV, and a two- 
component X, obtained on a 128 3 lattice, with the expan- 
sion of the universe assumed to be radiation dominated. 
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FIG. 1. Time dependence of the zero-momentum mode of 
<j> in units of its vacuum value for two runs with different 
realizations of random initial conditions for fluctuations. All 
other parameters were the same in both runs; r is conformal 
time in units of (A^ 2 (0))" 1/2 . 

Time dependence of the zero mode <^>o is shown in 
Fig. ||. Initially </>o oscillates with a large amplitude 
4> S> v. If all fluctuations were absent, the zero mode 4>q, 
in the expanding universe, would soon start oscillations 
near one of its vacuum values, ±i>. This would happen 
when the amplitude of the oscillations became smaller 
that V2v. In Fig. we see that the actual dynamics is 
completely different. The zero mode of the field </> con- 
tinues to oscillate near = even when its amplitude 
becomes much smaller than v. In other words, the field 
oscillates on top of the local maximum of the bare poten- 
tial. This can occur only because the effective potential 
acquires a minimum at (f> = due to interaction of the 
field with (X 2 ). 

At still later times, r > 720 in Fig. |, the zero mode 
of (f) decays completely. This should be interpreted as 
restoration of the symmetry <f> — > — tj> by nonthermal fluc- 
tuations. Finally, at r > 860 a phase transition occurs, 
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and the symmetry breaks down. 




transition. This is a signature of a bubble of the new 
phase and of the "soft" fluctuations produced in bubble 
collisions fllal. 
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FIG. 2. Variance of X (solid line), variance of cf> (dot- 
ted line), and amplitude of the zero-momentum mode of <f) 
(dashed line) as functions of time. 

Variances of the fields X and </> as functions of time 
are shown in Fig. ^. One can readily estimate that fluc- 
tuations of X (whose variance is multiplied by a larger 
coupling constant g 2 ) give a dominant contribution to 
the effective potential, compared to fluctuations of cj> it- 
self. At the moment of the phase transition (r « 860 in 
Fig. |2[), g 2 (X 2 ) is still larger than Aw 2 . This means that 
the effective potential still has a local minimum at 4> = 
even after the phase transition, which confirms that the 
transition is of the first order. 

To make sure that the phase transition itself is not 
a lattice artifact, we have varied the size L of the inte- 
gration box and the number of grid points. Values of L 
should be chosen in such a way that, on the one hand, 
there is enough control over the ultraviolet part of the 
spectrum for the results to be cutoff-independent and, 
on the other hand, the infrared part of the spectrum is 
represented well enough, so that the order of the phase 
transition is determined correctly. Good choices of L are 
made by monitoring power spectra of the fields, such as 
those shown in Fig. ||. 

At late times (but prior to the phase transition) the 
power spectra weakly depend on time and are power law 
functions at small k with an exponential tail at large k 
(conformal momentum k is defined in units of \/A</>(0), 
k = pa/ y/~\<fi(0)) . We made sure that the exponential tail 
of the power spectra is resolved in our simulations. Under 
that condition, the power laws do not significantly change 
with L, the number of grid points and, additionally, with 
the number of components of X. The power law for X 
is fitted by P x (k) cx k" 2 - 2 for k > I. At r w 700, 
k ss 10 corresponds to physical momentum p w gv . Later 
on, this value of the physical momentum corresponds to 
larger k, and at the moment of the phase transition it 
corresponds to momenta on the exponential tail of the 
power spectrum, which signifies that the condition (ii) 
is fulfilled, see Fig. || Note the enhancement of power 
spectra at small k occurring during and after the phase 
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FIG. 3. Power spectra of the fields X\/a (solid line) and 
4>/a (dotted line) before and after the phase transition. 

We have found that the transition grew stronger when 
we captured more of the infrared region, either through 
an increase of the grid size at fixed maximal momentum 
or through an increase of L itself. Thus we expect that 
the phase transition will be even more strongly first-order 
if we use larger lattices for our simulations. 

The transition also became more strongly first-order 
when we increased g 2 /X or the number N of fields Xi. 
Nevertheless, we observed a first-order transition even at 
N = 1. 

Direct demonstrations of the first-order nature of the 
transition are obtained by plotting a probability distri- 
bution of values of <f> over grid points and the actual field 
configuration in space. These are shown in Figs. ||, @. 
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FIG. 4. Probability distribution of <f>(x.) at several moments 
of time just before (dashed line) , during (solid line) , and at 
the end of the phase transition (dotted line). 

Prior to the phase transition, we observe a single peak 
in the probability distribution function of (f>. The peak is 
centered at cj) w 0, although its position oscillates slightly, 
synchronously with the oscillations of the zero mode. The 
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width of the peak is related to the variance ((<5</>) 2 ). At 
the moment of the phase transition, a distinctive second 
peak appears near the vacuum value of the field (position 
of this second peak is shifted somewhat from the vacuum 
value, because fluctuations still give a significant contri- 
bution to the effective potential). After that, the second 
peak grows, while the peak at <fi = decreases and even- 
tually disappears. During this period, positions of both 
peaks do not change. For some realizations of random 
initial conditions for fluctuations, the second peak ap- 
peared at <j) w v; for others, at <j> « —v. The time of the 
phase transition was also changing somewhat with the 
realization. This is an unambiguous signature of sponta- 
neous nucleation of a bubble of the new phase in a volume 
occupied by the old phase. 

By plotting the actual field configuration in space, we 
have indeed directly observed nucleation of a single bub- 
ble and the bubble's subsequent expansion until it was 
occupying the whole integration volume. The field config- 
uration at the beginning of the phase transition is shown 
m Fig. | To the best of our knowledge, this is the first 
time when lattice simulations allowed one to see nucle- 
ation of bubbles during a first-order phase transition. 



domains. Such models are not ruled out and in fact may 
have an interesting observable consequence, an enhanced 
(by bubble wall collisions) background of relic gravita- 
tional waves produced by the mechanism proposed in Rcf. 

. In models where the field </> has many components, 
the phase transition can lead to creation of strings or 
monopoles, instead of domain walls Jl5j ]. Finally, if the 
ratio g 2 /A is sufficiently large, one may expect a short 
secondary stage of inflation g] . To investigate this possi- 
bility one would need to study models with g 2 /X ^> 10 2 
or, equivalently, g 2 3> 10~ n , which is quite realistic. 
However, numerical investigation of this regime requires 
lattices of a much greater size than we currently use. We 
hope to return to this issue in a future publication. 
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AST95-29-225 (L.K. and A.L.), and by the Alfred P. 
Sloan Foundation (S.K.). 





FIG. 5. Bubble of the new phase. We plot a surface of 
the constant field 4> = —0.7v at the beginning of the phase 
transition. Inside the surface <j> < —0.7v. 

Let us describe some cosmological implications of our 
results. Models that exhibit behavior shown in Fig. |l| 
will lead to domain structure surviving till present. This 
conclusion is important because it allows us to rule out 
a large class of cosmological models that lead to domain 
wall creation. A different behavior of the zero mode is ob- 
served at smaller values of g 2 /A or N (but still g 2 /X 3> 1). 
There, the phase transition occurs when the zero mode 
still oscillates near <f> = with a relatively large ampli- 
tude. This is a new, specifically nonthermal, type of a 
phase transition. In such cases, bubbles of +v and -v 
phases will be nucleated in turn, but their abundances 
need not be equal, and for certain values of the param- 
eters one of the phases may happen not to form infinite 
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